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I. INTRODUCTION 

\Q ■ 

In nuclei having not-too-large neutron excess N — Z = 2T (T is the ground-state isospin value) the familiar giant 
fN) dipole resonance (GDR) is splitted into two components, characterized by the isospin values T< = T and T> = T + 1. 
q ' The GDR T>-component (GDR>) is the isobaric analog of the charge-exchange (in the /3 + -channel) giant dipole 

resonance (i.e. of the GDRM). Thus, the GDR> presents the specific double giant resonance (see, e.g. Ref. 
^| ' Because of Pauli blocking, the GDRW is not strongly collectivized and contains along with the main peak a few 
low-energy -states of the neutron- (proton-hole) type 0, These states is a base for the GDR> low-energy 
part, which consists of a few isolated dipole isobaric analog resonances (DIAR) 3]. One of these resonances in 90 Zr 
(E x = 16.28 MeV) has been identified in reactions with protons, electrons and photons (see Ref. [J] and references 
therein). The comprehensive study of the mentioned DIAR with the (e, e'p)-reaction has been reported in Ref. [H, 
where along with the form factor for two values of the transferred momentum the elastic El-radiative width and 
partial proton widths for decays into one-hole states of 89 Y have been deduced for this resonance. Being motivated 
by these experimental results, which are still not explained, in the present work we attempt to describe the gross and 
(direct) decay properties of isolated DIAR within a continuum-RPA(CRPA)-based approach. The basic elements of 
such description are presented in Sect. [TTJ The applications to the DIAR in 90 Zr are given in Sect. IIIH where the 
calculation results are compared with available experimental data. 
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II. CONTINUUM-RPA-BASED DESCRIPTION OF DIAR. 



o . 

, A. DIAR properties 

The dipole isobaric analog resonances in a target nucleus (Z, N) correspond to (2particle-2hole)-type excitations 
with a small admixture of (proton- (proton- hole) )-type states. Namely due to this admixture, El-decay of the DIAR 
into the target-nucleus ground state and also (direct) proton decays into one-hole states of the product nucleus 
(Z — 1, N) are possible. The description of the DIAR properties is simplified, provided that the isobaric analog state 
(IAS) is considered as the "ideal" one. The "ideal" IAS exhausts 100% of the Fermi strength and its wave function is 

generated by the Fermi operator — J2 a Ta . (The real IAR exhausts about 95% of the Fermi strength). In this 
approximation the wave function |s>, 1M) of a GDR> -component can be expressed via the parent GDR^ + ^-component 
wave function \s^ + \lM): 

|s>,lM) - (2T > )- 1 / 2 T ( -)| S (+),1A/). (1) 

In the mean-field approximation, the "ideal" IAS excitation energy equals the Coulomb displacement energy: = 
(2T>)~ X /[/ c (r)n(-'(r)(l 3 r, where U c (r) is the mean Coulomb field, vS \r) is the neutron-excess density of the 
parent nucleus (Z — 1, N + 1). (We assume that the IAS energy is independent of the parent state). 

Using Eq. (JTJ) , we express the GDR> strength function (corresponding to the single-particle external field V^\x) = 
(-1/2)t( 3 VYim) via the GDR( + ' strength function (corresponding to the external field V^ + \x) — t' + VYim): 

= (2T > )" I S(+)(4 +) ) 1 (2) 

where lu is the target- nucleus excitation energy, = lu — Aq 0]- The parent-nucleus dipole strength function 

5(+)((jj(+)) can be calculated within any CRPA-based approach (see e.g. Refs. @, Q). The cross section of photo- 
absorption, accompanied by GDR> excitation, is proportional to the strength function S>\u>): 
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where B = (167r 3 /3)(e /Tic). This cross section, integrated over an isolated DIAR (with the energy u> s ), (<7„ )*"*, deter- 



mines the El-radiative width for DIAR decay into the target-nucleus ground state: r 70 , s = (l/3ir 2 )(Lu s /hc) 2 (a> )* nt . 

Being applied to an isolated DIAR, Eq. |J2) allows one to express the (one-dimensional) particle-hole DIAR transi- 
tion density via the corresponding parent-state transition density: /3> )S (r) = (2T>)( -1 / 2 ) ps(r) . In particular, this 
transition density can be used for evaluation of the DIAR form factor, which determines the (ee')-reaction cross 
section [ij]. Within the impoved plane wave approximation the form factor equals: 



\F s {q)\ 2 = l2TiZ~ 2 y 2 {Z) 



P>,s{r)ji(qr)dr 
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where y(Z) = — — — ( 1 — exp ( — I I is the correcting factor (see, e.g., Ref. Q), and ji(qr) is the dipole 
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spherical Bessel function. 

The basic quantity in CRPA-based description of DIAR proton decay is the amplitude of the partial (7p)-reaction, 
accompanied by DIAR excitation. The poles in the energy dependence of this amplitude correspond to DIAR, while 
the pole residue is proportional to the product of the amplitudes of the elastic radiative and partial proton widths. 
Thus, the joint analysis of the corresponding calculated photo-absorption and (7p)-reaction cross sections allows one 
to evaluate the partial proton widths for a given DIAR. 



B. CRPA equations 

The CRPA expressions for the strength functions (uj^) and (7p)-reaction amplitudes can be obtained taking 
the CRPA equations in the form consistent with a phenomenological version of the Migdal's finite Fermi-system theory. 
As applied to description of charge- exchange giant resonances, these equations are given in detail in Refs. Q and 
for reader's convenience we use, in the main, the notations from these references. 

Of particular importance is the equation for the effective single-particle external field V^' (x, uj^). The dif- 
ference between the effective field and the corresponding external field V^(x) = V(r)YLM(n)T^' is due to the 
polarization effect from the particle-hole interaction taken for description of the non-spin-flip excitations in the form: 
F p - h -> 2F'(t(+M-)' + TH T (+)')i(f- 7). The effective-field radial part (r, oj^) determines all the main 
characteristics of a given charge-exchange giant resonance: the strength function Sy^u;^), the energy-dependent 
transition density r~ 2 p^(r, cj^^Ylmj and the nucleon-escape amplitude My](u>^) (c is the set of the nucleon- 
decay-channel quantum numbers). In particular, the dipole strength function S^(u>^ + ') of Eq. ^ can be directly 
calculated via the radial part of the corresponding effective field (see, e.g., Eq. (1) of Ref. [5(). Here we show the 
useful relationships, which have not been exploited in Refs. 0,(3, The pole representations of the effective-field 

radial part and strength function are valid in the energy region of each collective neutron- (proton-hole)-type state 
\s( + \ 1M) with the excitation energy below the neutron separation threshold: 

r - ojy' +i0 * u)M - ujV> + iO 

Here, p» (r) is the transition density of the (stationary) collective state \s^ + \lM), and di + ^ = J rpi + \r)dr is 
the dipole matrix element. In particular, this matrix element determines the partial radiative width Ty a of the 
corresponding DIAR fSubsect. Ill Ap . 

The transition-density elements, which define the transition density pi + \r) — P^n) s X-n{ r )Xv{f), can be found 
from the (formally, infinite range) system of the homogeneous equations, which is usually used within a discrete- RPA: 

PIV S = $ w Uu ~ n \ +) 2F' f P i + \r) Xn (r)xAr) %■ (5) 

Here, n v (n„) are the occupation numbers for neutrons (protons); r X v ( r Xw) are the radial bound-state wave 

functions for neutrons (protons); v — n rM ,j u ,l u is the set single-particle quantum numbers; \/3t[}J = ((v)\\Yi\\(tt)) is 

the reduced matrix element; [v) — j u , l v ; (ir) = j v , l v . Note, that the ratio —pm) s / (n„ — n,) = ci^ s is the amplitude 
of probability to find a certain 1~ neutron- (proton- hole) configuration in the collective state \s^ + \ 1M). 

The pole representation of the effective field, strength function, and also nucleon-escape amplitudes is valid in 
the energy region of an isolated doorway-state resonance corresponding to the respective (quasi-stationary) collective 
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state. As applied to the IAR with the excitation energy higher the proton separation threshold, the Breit-Wigner- 
type parametrization of the Fermi strength function Sp(u>^) and partial proton-escape amplitudes Mp y {uj^) are 
explicitly given by Eq. (|9|) of Ref. jfj|. (See also Eq. (6) of Ref. Q). We show here only the representation of the Fermi 
effective field: 

v ( C-h 2F ' S A 2 PA(r) ... 



Here pA(r) — y/4~xr 2 pa{t) is the IAS transition density; S]( 2 = y^n J pA(r)dr, Sa is the IAS Fermi strength; 
r = ^2 V T V is the IAR total proton width, while T v is the partial width for IAR proton decay into the neutron- 
hole state v^ 1 of the product nucleus. Note that the "ideal"-IAS transition density p % ^{r) = (N — Z)^ 1 / 2 rS~\r) is 
proportional to the neutron excess density. (S^ = (N — Z), the neutron excess is related to the parent nucleus). 



C. (7p)-reaction amplitude 

The expression for the amplitude of the (7p)-reaction, accompanied by DIAR excitation, is derived in three steps. 
First, we derive the CRPA expression for the amplitude of the IAR partial proton width, (r^) 1 / 2 . This expression 
can be obtained from Eqs. (8), (9) of Ref. @ (or from Eqs. (5), (6) of Ref. Q) and the pole representation of the 
Fermi effective field of Eq. © : 

{T v ) 1 / 2 = vViVy 2 5 (7r)W J x E (,)M<(^,(r)*, (7) 

Here, N v = (2j„ + l)n„ is the number of neutrons filling level v in a parent-nucleus state; r Xe(ir){ r ) is the normalized 
to the <5-function of energy radial continuum-state (real) wave function for the escaping proton with the energy 

e = e u + u>^; Vj[(f) — 2F'pA{r) is the IAS transition potential. For the "ideal" IAS w7' = Ac, and the transition 
potential (v^y 1 — (N—Zy' 2 v is proportional to the mean-field isovector part, provided that the isospin-selconsistency 
condition is fulfilled: v(r) — 2F'rS~\r), where v(r) is the symmetry potential (see, e.g., Ref. Q). 

The above-outlined CRPA expression of Eq. ([7]) can be applied to the IAS based on the "valence neutron + 
closed-shell core" parent-nucleus state. (In such a case N v = 1). This point allows us to generalize Eq. ([7]) for the 
case of the collective neutron- (proton- hole)-type parent states \s^ + \ 1M). Considering the IAS as the "ideal" one, 

we can substitute the product N^ 2 Xv(r) in Eq. (J7J) by the sum J2 V c tir Xf( r ) n ^ to get the expression for the 
doubly-partial amplitude of the DIAR proton width: 

(r^ ) J 1 / 2 = V2^(2T > )- 1 / 2n i/ 2 iW )w y' XtW (r)«(r) 9M (r/ )£ , + W W)2F'p(+)(r'^(r')^. (8) 

Here, e = + u>i + Ac is the energy of the escaping proton, the symmetry potential v(r) is defined for the 
(Z — 1, N + 1) parent nucleus, (rr')~ 1 g( 7r >)(r, r', e') is the radial neutron Green function with (y) — (it 1 ). In derivation 

of Eq. ([8|) we use the relation between the probability amplitude ci^ s and corresponding transition-density element 
of Eq. ([5|) and, also, the spectral expansion for the radial neutron Green function. The partial width for DIAR proton 
decay into the proton-hole state (7r) _1 of the product nucleus (Z — l,N) equals r* = Yl(n') ^(yw- Here, the sum 
is taken over the proton continuum-state quantum numbers (it 1 ) = which are compatible with the selection 

rules for the dipole transitions. It should be noted, that the use of Eq. ([7]) for derivation of Eq. ([8]) means, in fact, 
the use a perturbation theory on DIAR coupling to the single-proton continuum. 

Using the expression of Eq. |5]) for the doubly-partial DIAR proton-decay amplitude, one can obtain the CRPA- 
based expression for the corresponding doubly-partial (7p)-reaction amplitude R^n^w). For this purpose, we use 
the pole representation of Eq. ([4]) for the dipole effective field and dipole strength function and also the relation of 
Eq. ^ between the GDR> and GDR + strength functions. The result is: 

iV).M = (Bcj) 1 / 2 (2T > )- 1 n 1 JH^ t)M J x e ( ff oW«W^o(^ r ^ e - + w c +) )^(^^4 +) )x 7 ^(r')rfrdr^ (9) 

where e = £ n + u>. The pole representation of this amplitude is consistent with Eq. ©: 

fl ,,,,, M ^ ^>' /2 <W" (10) 

y 2 w s uj — uj s + iQ 
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Being based on Eqs. ([9]), (fT0|) , and the expression for the DIAR elastic radiative width fSubsect. HI A|) . we can evaluate 
within the CRPA the partial proton widths for a given DIAR. An expression similar to Eq. ([5J) has been derived in 
Ref. [|[ in an intuitive way. 



III. PROPERTIES OF THE DIAR IN 90 ZR 
A. Input quantities and model parameters 



A nuclear mean field and particle-hole interaction are the input quantities for any CRPA-based approach. Along 
with the phcnomcnological (non-spin-flip and momentum-independent) Landau-Migdal particle-hole interaction, we 
use a realistic phenomenological isoscalar part of the nuclear mean field (including the spin-orbit term), while the 
mean-field isovector part (the symmetry potential) and also the mean Coulomb field are calculated selfconsistently 
(see, e.g. Refs. ((Hsl). The parametrization of the mean-field isoscalar part is explicitly given in Ref. @. To describe 
better the neutron single-quasi-particle spectrum in 89 ' 91 Zr, the strength parameters are somewhat specified: 

U = 54.9MeV, U S o = 13.85MeVfm 2 , /' = 1.03, (11) 

where F' = f 300 MeV fm 3 . 

Following Ref. [H, we take effectively into account the contribution of the isovector momentum-dependent 
forces in formation of the GDR^) by the scaling transformation of the CRPA strength functions: S^(u>^') — > 
(1 + k')^ 1 (ut^) / (1 + k')). The parameter k! describes the relative contribution of the momentum-dependent 

forces to the corresponding energy- weighted sum rule. For 90 Zr the k! = 0.1 value has been found in Ref. [|| from 
independent data [10( . In the properly modified CRPA calculations of the photo-absorption cross section a> (oj) (see 
also Subsect. Ill A[) . performed in Ref. Q for 48 Ca and 90 Zr, a few isolated DIAR have been found. 

To improve description of the DIAR in 90 Zr, in the present work we take proton pairing approximately into account 
by simply substituting in CRPA equations the occupation numbers by the Bogolyubov coefficients v 2 . The same 
modification is used to calculate the neutron-excess density in order to keep the isospin-selfconsistency condition 
[HE!- The energy gap parameter is taken 1.0 MeV (this value is close to that used in Ref. [Hj]), while the proton 
chemical potential is properly calculated to find, finally, the coefficients v 2 . 



B. Calculation results 



We calculate first the photo-absorption cross section (oj) for the excitation energy interval 14-18 MeV, which 
corresponds to the low-energy part of the GDR( + ' strength function of Eq. © Q . In strength function calculations 
(performed with and without taking proton pairing into account) we add an artificial small imaginary part to the 
single-particle potential, used within the CRPA for calculations of the w-dependent quantities, to make presentation 
of the results more convenient (Fig. [1]) . Within the mentioned interval we find a few DIAR, whose energies uj s and 
calculated via (ct>)™* radiative widths r 70jS are given in Table HI The particle-hole part of the DIAR transition 

density p>, s (r) = (2T > )^ 1 / 2 ^ pi + \r) is next calculated via the residues at the poles in the energy dependence of the 
dipole effective field of Eq. The calculated transition densities are shown in Fig. [2] These densities are used to 
evaluate within the improved plane wave approximation the DIAR form factor of Eq. ([3]) . The results are shown in 
Fig. El Further, Eqs. ©, (|10p are used to evaluate the DIAR partial proton widths T^-^ for population of the 2piy 2 
(ttq), lgg/2 (ti), 2p 3 / 2 (7r 2 ) , I/5/2 fe) one-hole states in 89 Y. The calculated widths are given in Table U together 
with the experimental data of Ref. 0], related to the certain DIAR (E x = 16.1 MeV). For this DIAR we show also 
the calculated doubly-partial proton widths for decay into the main channels (ttq) and (7^) (Table ITT)) . 

Using the calculated DIAR total proton width T s = J^-n -T^ we can evaluate the (7Pt„t)-reaction cross section for 
the considered excitation-energy interval: 

= K>rr \ 2 (12) 



To compare with the experimental data of Ref. [ll|, we make "apparatus" averaging this cross section with the 
function I(ui, oj') = (I/tt)/((oj — oj') 2 +I 2 ). The averaged cross sections, obtained for the appropriate / values (/ = 50 
and 100 keV), are shown in Fig. 0] together with the experimental data of Ref. [l2j]- In calculations of the mentioned 
cross section the relatively small contribution of the GDR< excitation is also taken into account. It is done within 
the similar approach by the methods described in Ref. [lOj]. 
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C. Discussion of the results 

Due to the specific shell-model structure of 90 Zr, the proton pairing in this nucleus is rather weak. In particular, the 
number of protons, occupying the it level, N v = (2j Tr + 1) is relatively small for the l<?g/2 level. (In Table HT1 the N v 
values are shown for a few levels near the proton chemical potential) . Nevertheless, due to pairing, the DIAR with the 
energy 16.3 MeV appears in calculations (Fig. [I]) and DIAR proton decay into the 9/2 + -state of 89 Y becomes formally 
possible (Table [I]). Among the IAR found in calculations, three of them have relatively large elastic radiative width 
(Table H|, transition density (Fig. [2]) and, therefore, the form factor (Fig. [3]). The DIAR with the energy 16.1 MeV can 
be related to the resonance experimentally studied in Ref. [ij . Within the present approach we obtained a satisfactory 
description of the decay properties of this DIAR (Table JJ. The d 3 / 2 component of the proton partial widths is found 
to be the main one for decay into the 7rn-decay channel. This observation is in agreement with the analysis of the 
DIAR wave function performed in Ref. [4j. However, it is not truth for another main decay channel (7^) as it follows 
from Tables HI and HT1 The two-node transition densities of all DIAR (Fig. \2§ show that the low-energy 1~ parent 
states are not strongly collectivized, as it expected for the pygmy resonance. (The main- peak GDR(+) transition 
density is nodeless). The form factor, calculated for the DIAR with the energy 16.1 MeV, is in reasonable agreement 
with experimental data of Ref. Q . The calculation results show that at least two DIAR (with the calculated energy 
16.7 and 17.2 MeV) can be found experimentally via the (e, e'p)-reaction. This statement is partially supported by 
comparing the calculated and experimental (7p tot )-reaction cross section within the excitation-energy interval 14-18 
MeV (Fig. |4|). It is noteworthy that the approximations used in description of the IAS, proton pairing, and single- 
particle level scheme lead to uncertainties in evaluation of the DIAR energies (presumably, in a few hundreds keV). 
Description of the DIAR decay properties seems to be not much sensitive to these uncertainties. 

In conclusion, we first present continuum-RPA-based description of the decay properties of the specific doubly- 
collective states, 1~-IAR, and describe satisfactorily the available experimental data. 
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TABLE I: The energies (in MeV), elastic radiative and partial protons widths (in eV and keV, respectively), calculated for the 
isolated DIAR in 90 Zr. The parameters deduced in Ref. for the observed DIAR (the widths are given with errors) are also 
shown. 





"to 


r„- 






T-n-3 


14.91 


68 


49 


0.1 


8 




16.09 


110 


55.1 


1.7 


10.6 


0.12 


16.28 


108 ± 35 


50 ±6 


< 2.0 


20 ±5 


7± 2 


16.35 


66 


8 


6.5 


20 


0.10 


16.73 


169 


6 


0.4 


159 


0.02 


17.19 


185 


8 


0.3 


81 


0.06 



TABLE II: The doubly-partial widths (in keV), calculated for the main channels of proton decay of the DIAR at 16.1 MeV in 
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FIG. 1: The cross sections of photo-absorption, accompanied by excitation of the isolated DIAR in Zr. The full and dotted 
lines correspond, respectively, to calculations with and without taking the proton pairing into account. 
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FIG. 2: The calculated proton- (proton-hole) part of the DIAR transition density. The dashed, full, pointed, dashed-dotted and 
dashed-dot-dotted lines correspond to the DIAR at 14.9, 16.1, 16.3, 16.7, 17.2 MeV, respectively. 
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FIG. 3: The DIAR form factors, calculated within the improved plane wave approximation. The notations are the same, as in 
Fig. [2] The experimental data concerned with the DIAR at 16.1 MeV are taken from Ref. [fj]. 
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FIG. 4: The 90 Zr(7p tot )-reaction cross section, calculated for the excitation-energy interval 14-18 MeV (full line). The dashed 
and dotted correspond to averaging with the parameter I equal 50 and 100 keV, respectively. The thin line is for GDR< 
contribution to the cross section. The experimental data are taken from Ref. 



